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Abstract 

We use an analytic toy model of turbulent convection to show that most of 
the scaling regimes are spoiled by logarithmic corrections, in a way consistent 
with the most accurate experimental measurements available nowadays. This 
sets a need for the search of new measurable quantities which are less prone 
to dimensional theories. 

47.27 -i Turbulent flows, convection and heat transfer 
47.27.Eq Turbulence simulation and modeling 
47.27.Te Convection and heat transfer 

When an horizontal layer of fluid is heated from below, a heat exchange from the bot- 
tom to the top layer occurs. Dimensional consideration show that the non-dimensional heat 
exchange Nu depends on the geometry (via for example the aspect ratio), on the boundary 
conditions, and on the Rayleigh Ra and the Prandtl number Pr. Further dimensional anal- 
ysis of the dynamical equations governing the convective layer suggests that this dependence 
be a power-law. The exponent of the power law depends on the (dimensional) basic hypothe- 
sis. The classical theory predicts that Nu ~ Ra 13 with (3 = 0.33, which seems to be observed 
in electro-chemical convection [|IJ. Other theories lead to (3 = 0.2 or 0.24 [Q], (3 = 2/7 ||,|4]] 
or (3 = 1/2 || depending on the regime considered. A unifying synthesis of the possible 
scaling regimes, including their Prandtl dependence has been recently made by Grossman 
and Lohse ||[/]]. From the experimental point of view, the situation is rather unclear, since 
almost all the range of exponents between 0.25 and 0.33 has been measured (see || for a 
recent detailed review of the experimental measurements). Furthermore, recent experiments 
with fluids subject large variation of the Prandtl number have led to a new exploration of 
the phase parameter, and uncovered new scaling laws as a function of the Prandlt number. 
This new exploration can be used as a stringent tool to discriminate between the various 
theories, since two competing theories can provide the same value of f3, but different Prandtl 
number dependence. 

On the other hand, it is not quite clear whether the present dimensional theories sat- 
isfactorily capture the dependence of Nu as a function of Pr and Ra. In a recent serie of 
experiments, conducted with acetone at various aspect ratio, || show that their measure- 
ments are inconsistent with a single power law Nu(Ra), because the local effective exponent 
f3 e ff = dlnNu/dlnRa varies continuously with Ra in the range of Rayleigh number they 
consider. This effect can easily be accounted for in the scaling theory by considering a su- 
perposition of scaling laws, as proposed by . There is another solution, connected with the 
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existence of logarithmic corrections to scalings (see e.g. H). The difficulty associated with 
such a solution is that logarithmic corrections cannot be uncovered by dimensional analysis 
(nor numerical simulations!), and require detailed analytical computations. At the present 
time, they are no available analytical solutions of the full Boussinesq equations describing 
the dynamics of the convective motions, in the turbulent regime. Some time ago, Malkus 
H proposed to combine a weakly non-linear theory past the critical point, combined with 
a maximum principle to obtain approximate analytical solutions of the full problem. This 



approach developed further by Howard, Roberts, Stewartson and Herring (see refs |L0| for 
a summary) gives Nu = 0.24(i?a lni?a) 3//10 . Such a law with logarithmic correction fits the 
recent experimental data by Niemela et al. |TTJ very well. 

In the present work, we explore the predictions obtained using a solvable model of tur- 
bulent convection. This model couples large scale mean sheared velocity and temperature 
fields U = (U(z),0, 0) and Q(z), with small scale random velocity and temperature fields. 
This kind of large scale geometry is generally accepted as representative of the Boussinesq 
convection in the boundary layers, due to the shearing effect of the large convective cells. 
The model is closed by deriving an equation for the random component from the Boussinesq 
equation using two simplifying assumptions: i) the non-linear interactions of the small scale 
scales between themselves is modeled via a turbulent viscosity; ii) the small scale generation 
via the breaking of large scale structures is modeled by a random small scale forcing with 
prescribed statistics. This results in a linear stochastic equation for the random small scales, 
which can be analytically solved in the shear flow geometry by a decomposition of the small 
waves into localized wave-packets. In [0,0], the solution of the turbulent model were com- 
puted in the restricted case of 2D geometry, and turbulent Prandtl number equal to unity. 
This leads to expression of the mean and fluctuating velocity and temperature profile, as a 



function of Nu, Ra and Pr. In [113, these results were combined with dimensional estimates 



to derive the scaling laws governing the transport of heat. In the present contribution, we 
proceed one step further and use the analytical predictions about the profiles to compute 
directly the heat and energy dissipated across the turbulent boundary layer. From this, we 
derive an analytical expression of the non-dimensional heat flux, as a function of Pr and 
Ra. 

There can be a debate about whether such a simple model correctly accounts for the 
complex dynamics of the full, turbulent Boussinesq equations. A systematic numerical 
check of the various hypotheses pertaining the model has been undertaken in the past years 
TM and is still underway. In any case, we think that this model is very illustrative as to 



what kind of mechanism could take place in real, turbulent convection. Our computations 
show that most of the likely scaling regimes are spoiled by logarithmic corrections, in a 
way consistent with the most accurate experimental measurements available nowadays. We 
believe that this sets a need for the search of new measurable quantities which are less prone 
to dimensional theories. 

To set up notations, we start by a brief summary of the turbulent model [|l2"| , |l3f and of 
results obtained within it. In this model, the dynamics of the turbulent boundary layer is ob- 
tained from solutions of two coupled sets of equations. The first one described the dynamics 
of the mean (shear like) velocity U =< u >= (U(z), 0, 0) and temperature Q(z) =< 9 >: 

d t U + d z < uw > = -d x P + Prd 2 z U, 

d t e + d z < w'6' > = d 2 z e. (i) 
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These equations have been expressed, following Q in units the thermal diffusivity and the 
cell height. Here, Pr is the Prandtl number, the primes denote fluctuating (small scale) 
quantities and <> the averaging. In the stationary case, we get from ([!]) that d x P is a 
constant, independent of z. In the laminar case where < u'w' > and < w'& > are negligible, 
we thus obtain the well known parabolic profile for the velocity and the linear profile for the 
temperature. In the turbulent case, the profiles are linear within the thermal or the viscous 
layer, while outside this layer, they are given by the condition 

d z < u'w' >= -d x P, d z < w'6' >= 0. (2) 

To close the system, we need < u'w' > and < w'6' >. They are obtained as solution of a 
linear, stochastic equation valid for localized wave-packets of velocity and temperature: 

D tUi = -ihp - wd z U5 a + RaPr65 i3 - iVk 2 ^ + ft u) 

D t e = -wd z e-k 2 e + f ie \ (3) 

where 

u(x, k, t) = J g(\x - x'|)e ik -( x - x ')u(x / , t)<bd, (4) 

g being a function which decreases rapidly at infinity. We have dropped primes on fluctuating 
quantities for convenient notations and introduced the total derivative D t — df + Ud x — 
d z (Uk x )dk z ■ Note that the linear part of (|3]) is exact and describes the non-local interactions 
between the mean and the fluctuating part. The major approximation of the model is to 
lump the non-linear terms describing local interactions between fluctuations into a turbulent 
viscosity, or equivalently, into a turbulent Prandtl number Pr 1 . The forces /w and 
appearing in (^) are small scale random forces which are introduced to model the seeding 
of small scales by energy cascades (for example via turbulent plumes, detaching from the 
wall). 

The analytical solutions of ([]]) and (0) have been obtained in |r2|JT3|| in the 2D case 



(no movement in the direction transverse to the mean flow), and for Pr 1 = 1. There are 
numerical and analytical evidence that 2D geometry is sufficient to capture the physical 
mechanism responsible for the Nu scaling with Ra [ITS] and also to capture the correct 
shape of the equilibrium profile in the neutral case [0. Note that the vortex stretching, 
which is theoretically absent in 2D geometry, has been implicitly accounted for via the 
turbulent viscosity. The unit value of Pr* was dictated by our ignorance of the exact value 
of this parameter (which is likely to depend on Pr). We shall keep in mind that it may 
induce a wrong dependence of the correlations with respect to the Prandtl number. In any 
case, we shall consider these solutions as a toy model of turbulent thermal convection, and 



investigate their scaling properties. For the mean flows, the solution are |T3 



d z U ~ d z Q ~ \, (5) 



z z z 



resulting in a constant (with z) Richardson number: 



[d z U) 
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In developed convective turbulence, this number is large and negative Ri <C —1. This 
motivates a large Ri expansion of the solution of (||), which provides the expression of the 
correlations as : 



<wu > ps —— oc z, 
d z U 

< W 'Q> > ~ <^E^L oc cte = Nu 

RaPr 

ul 



<w' 2 > 



d z U' 



nil n U l9 z U uld z Q 

< > p^ — = —— ( 7 1 

\RaPr) 2 RaPrd z U' {> 

where we have introduced ul = [d z U /k 2 ) 2 ^~ Ri ^ . Note that the last equation of ([?]) just 
reflects the balance between the vertical energy and the buoyancy force. Such a balance can 
be expected to hold only at Pr > 1 0. For low Pr, temperature fluctuations were found to 



behave like a passive scalar, except in the boundary layers, or in the center of the cell |L7 
Since our expression pertain only to the boundary layers, we feel safe to consider it for both 
the Pr > 1 and the Pr < 1 regime. 

The second equation of (0) can be used to express w* as a function of Nu and Ri. Then, 
we can eliminate Ri by matching the turbulent profiles with the viscous or diffusive solutions 
PrdJJ = u T and d z Q = Nu, where u T is the friction velocity and Nu the Nusselt number. 
For this, we introduce the thermal length scale Xg = l/2Nu and the viscous length scale 
X u = KPr/u T (K ~ 0.4 is the Karman constant), and operate the matching via smooth 
functions according to: 

ulPr-- 1 

aJJ 
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The shape of d z U is dictated by the exact analytical solutions of the profile found by [TB 
in a neutral layer. The shape of d z Q was chosen as the simplest smooth function matching 
the viscous and turbulent layers. Using these expressions, we can eliminate Ri and write for 
the fluctuations: 



< w ' 2 > = Nu 3/2 (RaPry 



/2 



U- 



Nu 3/2 (RaPrf 12 — Jl + (zl\ u f 



< e' 2 > = u 



Nu 3/2 d,Q 



T {RaPrf/ 2 d z U 
Pr Nu 5 / 2 \fi + (z/X u y 



(9) 



u T (RaPr) 1 ' 2 1 + (z/X e ) 2 ' 

Note that these fluctuations exist only for z > X u . The expression we derived are only 
valid inside the turbulent layers, which we define as the locations were the fluctuations of 
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horizontal velocity and of temperature reach a maximum. The scales at which these maxima 
are reached are labeled Xblv an d Xblt, respectively. Beyond these scales, the flow enters 
the "bulk region" , in which the velocity and temperature fields are organized into convective 
large scale, whose details are non-universal and strongly depend on the geometry (aspect 
ratio) of the cell. It is difficult to evaluate analytically the vertical dependence of the profiles. 
On the other hand, support from data can only be provided by numerical simulations, 
encompassing relatively low Rayleigh numbers. One important outcome of these simulations 
p"8|^20| is the existence of two distinct flow regimes: one for Pr < 0.35 in which the flow is 
dominated by a large re-circulation cell, and one Pr > 0.35 in which isolated thermal plumes 
can develop. The detailed scaling of the velocity pattern has only been partially studied 
by Kerr |2lJ in the regime Pr = 0.7, Ra = 8 x 10 7 and for aspect ratio 1. He found no 
indication of the logarithmic region, but rather a variation of the horizontal velocity profile 
consistent with a z" 1 ^ 2 law. This can be easily understood by noting first, that in this Prandtl 
regime, the Reynolds number corresponding to Ra = 8 x 10 7 is too low to allow for a well 
developed turbulent boundary layer. Therefore, the logarithmic region is non-detectable. 
On the other hand, in the core, there is a tendency for strongly homogeneous temperature 
(due perhaps to the presence of plumes, which favor the good mixing), with nearly constant 
value Oo throughout the central region. The classical free fall velocity argument then predicts 
vertical velocities varying like W 2 ~ RaPrQ Q z, resulting, by incompressibility d z U ~ d z W 
in a z~ l l 2 law for the horizontal profile. Given the small extent of the logarithmic zone in 
this regime, we shall therefore approximate the bulk velocity in this region, at Pr > 0.35 by 
u T (z I 'A-u) -1 / 2 (to ensure the matching with the viscous layer). Unfortunately, no comparable 
measurements are available at low Pr. On the one hand, we can expect the logarithmic 
region to be more extended in this regime, because of the larger Reynolds number involved. 
Interestingly enough, inversion of the temperature profiles in the central region have been 



detected in the experiment of ||17 |, with a nearly constant and positive value of dT/dz 
(whether this is correlated with the absence of plumes is an interesting open question). In 
such a case, the free-fall argument (which is valid provided we are in the Bolgiano regime, i.e. 
at large enough scales) predicts vertical velocities varying like W 2 ~ RaPrz 2 , i.e., via the 
incompressibility, constant horizontal velocity profiles. Flatter profiles of horizontal velocity 
have indeed been observed in some numerical simulations of Kerr but the scaling has 
not been properly checked. In any case, as a first approximation, we shall assume that the 
horizontal velocity in the central region and in the Pr < 0.35 regime, is constant Uq. The 
matching of this constant with the logarithmic profile then imposes Uq = Ku T (\n(X ) + B), 
where A Q is the outer scale, and B is a constant which may depend on the Prandtk number 
(under neutral condition, B « 5). In the sequel, we assume U = fu T , and assume that / is 
independent of the Prandtl number (for Pr < 0.35). 

These estimates can be further used to determine the scaling behaviors of the two scales 
Xblt and Xblv characterizing the turbulent boundary layer. For this, we follow the logic of 
our turbulent model which states that in the boundary layer, the fluctuations are passively 
advected by the large scale horizontal velocity and obeys Ud x u' = Prd 2 u' or Ud x 8' = d 2 9' . 
On dimensional ground, this shows that the two scales are determined via the two implicit 
equations: 

1 Pr 
U(Xblt) ~ to— , U(Xblv) ~ Ta— - (10) 

"BLT ^BLV 



5 



Here, we have ignored the aspect ratio dependence. Aspect ratio obviously affects the scaling 
regimes, but since we cannot quantify its effect in the large scale circulation, we shall ignore 
it hereafter. Depending on the Reynolds number Re = u T /Pr and on the Prandtl number, 
different regimes of velocity will be probed by this implicit relations, and the scaling of Xblv 
and Xblt as a function of Pr and u T will vary. Introducing e(A) = d\nU/dhiz\ z= \ (where 
A = Xblv or Xblt), we can summarize this dependence as: 




1 x S(e)/(2+e) 



Xblt = Pr^+^u^^ 

V m Xblt 

Here, 5(e) is equal to 1 if e = and zero elsewhere. The dependence of e in Xblt and Xblv 
has been omitted for simpler notations. Also, the logarithmic and constant velocity regime 
have been lumped into the single notation U(X) ~ w T /(l/ In A + 1//) which patches the two 
consecutive characteristic behaviors. 

These considerations can be used to compute the Rayleigh and Prandtl dependence of 
the Nusselt number. For this, we follow and use the two rigorous relations involving 
the global average of the kinetic and thermal dissipation rates: 

< PridiUj) 2 > v = RaPr(Nu - 1), 

< (di6) 2 > v = Nu. (12) 

In ( |i2|) averages are taken over spatial volumes. Dissipation takes place both in the boundary 
layers, and in the "bulk" region. However, several independent observations militates in favor 
of the importance of boundary layers, rather than bulk region, to determine the scaling of 



the Nusselt with the Rayleigh number. For example, Julien et al [22| observe that when 
switching from no-slip to stress-free boundary conditions, the exponent of the power-law 
Nusselt versus Rayleigh changes from close to 2/7 to close to 1/3. More recently, Ciliberto 
and Laroche |23j observed that both the prefactor and the exponent of the power-law changes 
when switching from smooth to rough bottom plate in the convective cell. On the other 



hand, blocking the large scale circulation does not affect the heat transport [24]. In numerical 



simulations at Pr = 0.7 and Ra = 8 x 10 7 , Kerr [J2T] noted that a change of aspect ratio 



(which changes the large scale circulation) only affects the energy dissipation in the bulk of 
the flow. In the boundary layers, the total dissipation seems to remain always about 1/4 
of the total dissipation. Gathering all these information, we shall assume that a constant 
fraction of heat and energy is dissipated in the boundary layers, independently of geometry, 
Rayleigh number, Prandtl number... This allow us to consider similar equations to (ff2j), 
where the average is only taken over the volume encompassed by the boundary layers. The 
advantage is that, in that region, we can use our exact analytical solutions to compute the 
kinetic and thermal dissipation. 

In each case, these dissipation include two contributions: one from the mean flow, and 
one from the fluctuations. For the mean flow contribution, we use (p)) and find: 

3 , / Xblt \ 



e v = Pr f BLV dz(d z U) 2 = Ku 3 T arctg 
Jo 



X u 
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dz{d z Qf 



Nu 2 
~2X~o 



XbLt/X^ 



1 + (XBLr/Xey 



+ arctg 



'X 



BLT 



A 6 



(13) 



For the fluctuations, we first note that the velocity and thermal fluctuations are non- 
negligible only for z > X u Q Further, we use incompressibility to get (d z u) 2 ~ X B ^ LV {d x u) 2 ~ 
Xj^ LV (d z w) 2 . Then, we approximate d z w ~ d z \J < w 2 > and d z 9 ~ d z \J < 9 2 >. Using @ 
and performing the integration, we finally obtain for the contribution of the fluctuations: 



(NuPrf^Ra 1 / 2 



{Xblv/ A u ) 



AK 



(XbLV 

„2 



IK? 

2?/ 



+ In 



2JI + {Xblv / X u ) 2 + 2(A BLy /A u ) 



2v/2 + 2 



rfx 2 ) 2 



Vb£t/A u 

4K{RaPry/ 2 Ji (1 + x 2 f/ 2 {\ + r/ 2 x 2 ) 3 : 



f] = A„/A 6 



(14) 



Based on these expressions, we can singularize six different regimes. In the first three regimes, 
the dissipation is dominated by the mean flow contribution. This regimes are expected to 
hold at low Reynolds number (the turbulence is weak) or when the Prandtl number is low, 
according to the numerical study of |[L8|| . The thermal dissipation then fixes Xblt ~ Nu^ 1 . 
Then we can distinguish two cases: at very low Reynolds number, the boundaries are non- 
turbulent, and we have Xblv < ^u, so that e v ~ u\Xblv i 1 P r "> an d the two boundary layer 
scales follow ([TT]) with e = 1. In this regime, we get: 



Nu ~ Ra^Pr- 1 ' 12 , Re ~ i?a 3 / 8 Pr" 5 / 8 . 



(15) 



This Nusselt versus Rayleigh dependence corresponds to regime I u of 0. They argue that 
for fixed Raleigh number, it could explain the weak decrease of the Nusselt number observed 
at large Prandtl number. For larger Reynolds number, when Xblt > X u , we get to a regime 
where ejj ~ ~ RaPrNu. This corresponds to a Reynolds number varying like (RaNu) 1 ^ 3 , 
a situation observed in the low Prandtl number experiment of @. The final expression of 
Nu depends on how the thermal boundary layer scale Xblt ~ Nu^ 1 varies with the mean 
flow. Using (|TT|) , we obtain: 



Nu 



^ a (l+e)/(5+2e) p r (l-2e)/(5+2e) 
(l/ln(Mx) + l//)3*( £ )/(5+2 £ ) • 



(16) 



For e ~ 1 (which may be obtained at low Rayleigh, for Pr > 0.35), this corresponds to 
the 2/7 law of in which the Nusselt decreases like Pr^ 1 ^ 7 and the Reynolds number 
varies like Ra~ 3 ^ 7 . For low Prandtl, we have to consider the case e = 0. In that case, we 
observe an interesting regime, where the Nusselt varies like Ra 1 ^ 5 '(In Ra) 3/ ' 5 at moderate Ra, 
then turning into a Ra 1 / 5 regime for larger Ra. The logarithmic correction increases the 
scaling exponent up to a value close to 1/4, as shown in Fig. 1. In this regime, experimental 
measurements will tend to show first a 1/4 regime, followed by transition to a 1/5 regime, 



Temperature fluctuations are generated through velocity fluctuations, and, thus, are slaved to 
them 
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when the logarithmic corrections vanish, i.e. when the scaling of the boundary layer scale is 
determined by the bulk flow. This regime therefore exactly explain the experiments by ||. 

Other interesting regimes appear when the fluctuations dominate the kinetic or thermal 
dissipation. The conditions under which this situation appears cannot be determined by a 
mere estimate of the ratio of e u versus ejj, for example, because there are some indications 
that this condition is in fact subject to boundary conditions. For example, experiments 
conducted at same aspect ratio, same cylindrical geometry and same Prandtl number by 
|| and [^5| lead to different values of the Nusselt number. In ||, the scaling relation 
Re ~ (RaNu) 1 / 3 (characterizing the regime in which kinetic energy dissipation is dominated 
by the mean flow ) is satisfied. In the experiment of [P5|1 , a careful study of the length 
scales revealed a significant deviation from the law Xblt ~ 1/Nu which would be obtained 
in the regime explaining the experiment by [0, and in which the thermal dissipation is 
provided by the mean temperature profile. It is then logical to attribute the difference of 
Nusselt via a difference in boundary conditions, which would favor or not the dominance of 
fluctuations. For example, the thickness or composition of the bottom plate may play such 
a role (Castaing, private communication). 

Let us consider first the regime in which mean flow dominates the thermal dissipation, but 
fluctuations dominates the kinetic dissipation. In such a case, we still have Xblt ~ 1/Nu, 
but the scaling of Nu depends on the behavior of Xblv- Assuming (|TT|), and taking e = 1 
for Xblt (no fluctuations of temperature), we get: 

„ Ra^l^Pr^) [ln(AWA,)r 2(2+£)/(8+7€) n?1 
(l/ln(A Bi ,) + l//)-«(W • [U) 



For e = — 1/2 (valid for Pr > 0.35), we obtain two different behaviour with (17): one, with 
Xblv/Xu ~ 1- This situation prevails for example under stress-free boundary conditions, 
like studied by WernePC, or for very large aspect ratio. It may also occur at very large 
Prandtl number, if, as discussed in || the viscous length scale approaches the size of the 
box and cannot increase any further. In such situation, we get Nu = Ra^Pr^ 9 . In 
the opposite situation, where Xblv/X u 3> 1, we then have (within logarithmic corrections) 
Xblv ~ Re 1/3 . This gives Nu ~ Ra 1 / 3 / (\n(Ra)) 2 / 3 . This law reproduces extremely well the 



recent experimental measurement of Niemela et al. [PH , obtained with Helium (Pr = 0.7) 
over 11 decades of Rayleigh numbers. These measurements can be fitted by a power-law with 
exponent 0.308 over this range of Rayleigh number. However, as shown in Figure 2., even 
over such a wide range, the law with logarithmic corrections cannot be distinguished from the 
power-law exponent. Other measurements are also explained by this regime, for Pr > 0.35. 
For example, in this regime, we predict a variation of the Nusselt number with the Prandtl 
numer going like Nu ~ Pr -1 / 9 / m (Pr) (because Re depends on the Prandtl number). This 
variation is shown in Figure. 3. A power-law fit of this variation over 0.7 < Pr < 100 gives 
an exponent —0.2, exactly like in the experiment of Ashkenazi and Steinberg ||26|| . Also, the 
Reynolds number in this regime varies like Re = J Nu 3 / Pr. The power-law dependence 
mimicked by the logarithmic corrections is (see Fig. 4 and 5) Re ~ Ra 0A6 Pr~ 08 , again, 



in excellent agreement with the measurement of pq| . Other quantities, like the scaling 
of the central fluctuating temperature, or the length scales (Table IV) are also in excellent 
agreement with direct measurements of Niemela et al. [|ll|] or by Kerr Jl9| . We do not discuss 



the regime corresponding to e = in (fl7|) because, apparently, it is not applicable to the low 
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Prandtl regime: it assumes that Xblt < A u , at variances with experimental observations by 



Segawa et al |17 



The last regime is obtained when temperature fluctuations dominate the thermal dissi- 
pation, making Xblt differ from 1/Nu. The two interesting regimes then depend on whether 
r] = \ u /\e = Nu/Re is less or larger than 1. In the first case, we get: 

Nu 5 ^ 2 P^blt/K x 2 



4K(RaPry/ 2 Ji (1 + x 2 ) 3 / 2 

Nu 5/2 

~ 4K(RaPryrM ^r/X,l). (18) 

We then get from (0) that Nu ~ ( J RaPr) 1 / 3 /(ln(A BLT /A u )) 2 / 3 . This regime is similar to 
regime (|TT|), with a different prefactor. In the case where r] 3> 1, we get 

Nil 5 / 2 f^BLT/K dx 



4K(RaPry/ 2 Ji r] 2 x 2 
Nu^ul 



(19) 



4K(RaPry/ 2 ' 

From (|I~2"D, we get Re = u T /Pr = Pr~ 3 ^ 4: (RaNu) 1 ^. Using the expression (pTTJ ) for \blv-, 



and e u = RaPrNu, we finally get : 

^1/(3+2^(1+20/(3+20 ln(ABLv/Au) - 2 (2+ e )/(3+2 e ) 

(1/ H*BLV) + l//)-^)/(3+ 2 ^) • (20) 

For e = (corresponding to low Prandtl number), we obtain a regime Nu = Ra 1/3 with 
logarithmic corrections which depend on the Rayleigh number: for low Rayleigh num- 
ber, the scaling of the boundary layer is determined by the logarithmic region, making 
Nu ~ Ra 1 ^ 3 / ln(i?a) 8 / 3 . This regime mimics a scaling exponent close to 1/4 (see Figure 
6). For larger Rayleigh numbers, the scaling of the boundary layer scale is determined 
by the bulk flow, changing the correction to Nu ~ Ra 1 ^ 3 / (ln(Ra)) 4: ^ 3 . This results in a 
larger effective exponent, of the order of 0.28. In this regime, the boundary layer scale like 
Ablv ~ {Re\n.{Re))~ 1 / 2 , providing an approximate power-law in Rayleigh with exponent 
—0.2. All these features were observed in the mercury experiment of J27|. For e = —1/2 



(Pr > 0.35, we obtain a regime in which the Nusselt depends only weakly on the Prandtl 
number Nu ~ Ra 1 ' 2 /(ln(i?e)) 3 / 2 . This regime corresponds to the "ultra-hard" convective 
regime predicted by Kraichnan ||, but with logarithmic corrections. In Fig. 7, we compare 
this regime with the Nu ~ p r o.o72^> a o.389 approximate power-law measured by Chavanne et 
al [28| in Helium, at Ra > 10 11 . Our formula predicts a very weak dependence of in the 



Prandtl number, like in |28|, but with opposite sign. This could be accounted for minute 
variations of the mean profile around the value e = —1/2. However, the plus sign of the fit by 
Chavanne et al. could be an artifact due to the scatter of the data. Indeed several measure- 
ments performed at fixed Rayleigh number and for increasing Prandtl show that Nu tends 
to decrease rather than increase, as in the formula by Chavanne et al. The predictions in this 
regime are also in agreement with their Reynolds number measurements Re ~ Ra 0A9 Pr~ 75 
if one considers that this Reynolds number is based on the vertical velocity, i.e. Re ~ w c Pr. 
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Interestingly enough, the slight discrepancy between our prediction for d In w c /d In Pr (little 
less than 0.25) and the Chavanne et al. fit (little more than 0.25) can be traced back to 
the discrepancy in the Nusselt number. So there must be some systematic error present at 
this level, either in the theory, or in the experiments. In any case, this new regime may 
explain to puzzling facts: 1) why Glazier et al do not observe transition towards a new 
scaling regime, despite their very large Reynolds number. This is because they were already 
in the ultra-hard regime! 2) why for the same geometry and the same aspect ratio, Niemela 
et al did not observe a transition towards the regime observed by Chavanne et al. It is 
because, somehow, they do not allow the temperature fluctuations to grow enough so that 
their contribution supersedes that of the mean flow. 

A good summary of our results can be found in the tables. For completeness, we have 
indicated the scaling behavior of various quantities measured in experiments, like A c or w c 
(the value of temperature or velocity fluctuations at the center of the cell). These values 
can be estimated from ([]), with z = 1. When available, we have indicated comparison 
with measurements. This summary shows that many scaling behaviors observed in recent 
experiments can be reproduced using log-corrected scaling laws. Overall, the agreement 
of the theory based our toy model of thermal convection is excellent for the case Pr ~ 1 
(Helium). It predicts reasonably well the Rayleigh dependence of the quantities for low 
Prandtl number, but it is not clear whether the Prandtl number dependence in these regimes 
is well predicted. This might be caused by our initial hypothesis (taking Pr 1 = 1). Also, 
we were unsure where to put in this picture the recent experiments by || performed in 
water. The difficulty is that these experiments are performed at many different aspect ratio, 
and it is difficult to disentangle the aspect ratio dependence of the local scaling exponents 
with other dependences. At least in the case of Helium, apparently, we reached a point 
where experimental measurements cannot discriminate between logarithmic corrections to 
scaling, or corrections due to superposition of power laws, like proposed by ||. This calls 
for identification of more stringent tests of theories. One possibility, which has been so 
far under-exploited, would be to concentrate on statistical properties rather than average, 
i.e. on probability distributions functions. In this context, it is interesting to note that 
our toy model of turbulent convection allows such investigations, via a Langevin formalism 
exploiting the linearity and the stochasticity of the equations. In neutral condition, our 
toy model leads to PDF's for velocity with qualitatively similar behavior than the PDF's 
measured in high Reynolds number flows |29[]. The generalization of these results to the 
convective regime is underway. 

I thank the Fluid Mechanics group of the ENS Lyon for many very interesting discussions. 
My special thanks goes to B. Castaing for his continuous interest and support in the present 
work. 
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TABLES 



Q(X,Y) 


This theory 


Experiments 




dlnQ/dX 


dlnQ/dY 


dlnQ/dX 


dlnQ/dY 


Nu 


1/4 


-1/12 






Re = u T /Pr 


3/8 


-5/8 






A c 


-3/16 


-13/48 






w c 


7/16 


3/16 








-1/4 


1/12 






^BLV 


-1/4 


5/12 







TABLE I. Summary of local exponents of different physical quantities in regime 1: low 
Reynolds number, mean flow dominates. In this table X = In Ra and Y = InPr. No dependence 
on X or Y of the scaling exponent indicates real scaling with respect to Ra or Pr respectively. 



Q(X,Y) 


This theory 


Experiments 




dlnQ/dX 


dlnQ/dY 


dlnQ/dX 


dlnQ/dY 


Nu 


2/7 


-1/7 


0.282 




Re = u T /Pr 


3/7 


-5/7 


0.43 




A c 


-5/28 


-2/7 


-0.147 




Wc 


13/28 


1/7 


0.491 






2/7 


-1/7 








np 


np 







TABLE II. Summary of local exponents of different physical quantities in regime 2: Pr > 0.35, 
mean flow dominates. This regime corresponds to the theory of Q|. The measurements are those of 
the Chicago group ||. However, this experiments might well be better described by regime 4, see 
below. In this table X = In Ra and Y = In Pr. No dependence on X or Y of the scaling exponent 
indicates real scaling with respect to Ra or Pr respectively. The symbol np means not predicted 
by this theory. 
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This theory 


Experiments 




a in C^/ aA 


a m t// a i 


a In LI 1 aA 


A 1-^ /^7V^ 

a In / a i 


Nu 


l/5 + 3//5(/*+* 2 ) 


i/5 + 3//5(/y + y 2 ) 


0.26 then 0.2 


0.21 


Re = u T /Pr 


2/5 + //5(/* + X 2 ) 


-3/5 + //5(/y + y 2 ) 


0.424 




A c 


-l/5 + 3//20(/X + X 2 ) 


-i/5 + 3//20(/y + y 2 ) 






w c 


2/5 + 9//20(./* + * 2 ) 


2/5 + 9//20(/Y + Y 2 ) 








-l/5-3//5(/* + * 2 ) 


-l/5-3//5(/Y + Y 2 ) 








rap 


rap 







TABLE III. Summary of local exponents of different physical quantities in regime 3: Pr < 0.35, 
mean flow dominates. The measurements are from Q. In this table X = In Ra and Y = In Pr and 
/ is a constant, depending on the shape of the bulk velocity profile. This constant is not predicted 
by the theory. No dependence on * or Y of the scaling exponent indicates real scaling with respect 
to Ra or Pr respectively. The symbol np means not predicted by this theory. 



Q(X,Y) 


This theory 


Experiments 




dhxQ/dX 


dlnQ/dY 


dlnQ/dX 


dlnQ/dY 


Nu 


1/3 - 2/(3*) 


-1/9 - 2/3Y 


0.309 1 


-0.2 2 


Re = u T /Pr 


1/2 - 1/X 


-2/3 -l/Y 


0.43 2 ' 3 


-0.75 2 


A c 


-1/6-1/6* 


-5/18 - 1/6Y 


-0.145 1 




w c 


1/2 - 1/2A 


1/6 - 1/2Y 


0.53 4 


0.08 4 




-1/3 + 2/3* 


1/9 + 2/3Y 


-0.28 3 




^BLV 


-1/6 + 1/3* 


1/9 + 1/3Y 


-0.14 3 


0.21 3 



TABLE IV. Summary of local exponents of different physical quantities in regime 4: Pr > 0.35, 
fluctuations dominates. The measurements are from M (superscript 1), |2(| (superscript 2), [20[| 
(superscript 3) and (superscript 4). In this table * = lni?a and Y = InPr. No dependence 
on * or Y of the scaling exponent indicates real scaling with respect to Ra or Pr respectively. 
The symbol np means not predicted by this theory. 
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Q(X,Y) 



This theory 



dlnQ/dX 



dhxQ/dY 



Experiments 
dhiQ/dX klnQ/dY 



Re 



Nu 

= u T /Pr 
A c 

w c 

XbLT 
XbLV 



1/3 - 4/3X - 4//3(/X + X 2 ) 
1/3- 1/3X - f/3{fX + X 2 ) 

-1/Q-1/3X - f/3(fX + X 2 ) 
l/2-l/3X-f/3(fX + X 2 ) 
np 

-1/6 + 1/6X + f/6(fX + X 2 ) 



i/3-4/(3y)-4//3(/y + y 2 ; 

-2/3 -1/3Y - f/3(fY + Y 2 ) 

-i/6-i/3y-//3(/y + y 2 ) 
i/ 2 _i/ 3 y-//3(/y + y 2 ) 

np 

-i/6 + i/6y + //6(/y + y 2 ) 



0.25 then 0.285 



-0.2 



TABLE V. Summary of local exponents of different physical quantities in regime 5: Pr < 0.35, 
fluctuations dominates and Xg S> A n . The measurements from [27]. In this table X = \nRa and 
y = InPr and / is a constant, depending on the shape of the bulk velocity profile. This constant 



is not predicted by the theory. No dependence on X or Y of the scaling exponent indicates real 
scaling with respect to Ra or Pr respectively. The symbol np means not predicted by this theory. 



Q(X,Y) 


This theory 


Experiments 




dlnQ/dX 


dlnQ/dY 


dlnQ/dX 


d In Q/dY 


Nu 


1/2 - 3/(2X) 


-3/2Y 


0.389 


0.072 


Re = u T /Pr 


3/8 - 3/8X 


-3/4 - 3/8Y 






A c 


-1/8- 3/(8X) 


-1/4-3/8Y 






w c 


5/8- 9/8X 


1/4-9/8Y 


0.49 


0.28 


Xblt 


np 


np 






Xblv 


-1/8 + 1/8X 


1/4+ 1/4Y 







TABLE VI. Summary of local exponents of different physical quantities in regime 6: Pr > 0.35, 
fluctuations dominates and Xg S> X u . The measurements are from []28| . In this table X = In Ra 
and y = In Pr. No dependence on X or Y of the scaling exponent indicates real scaling with 
respect to Ra or Pr respectively. The symbol np means not predicted by this theory. 
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FIGURES 

FIG. 1. Nusselt vs Rayleigh in regime 3 (Pr < 0.35, mean flow dominates). The symbols are 
the power-law fits of the experimental measurements by Q. The line is the theoretical formula 
predicted by the toy model (eq. (|l6|) with e = 0). 
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FIG. 2. Nusselt vs Rayleigh in regime 4 (Pr > 0.35, velocity fluctuations dominate but tem- 
perature fluctuations are negligible). The symbols are the power-law fit Nu ~ i?a ' 309 of the 
experimental measurements by The line is the theoretical formula predicted by the toy model 
(eq. (0) with e = -1/2). 
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FIG. 3. Nusselt vs Prandtl in regime 4 (Pr > 0.35, velocity fluctuations dominate but temper- 
ature fluctuations are negligible). The symbols is the theoretical formula (eq. (|l~7|) with e = —1/2). 
The dotted line is a power-law fit, with slope —0.2, mimicking the experimental fit of [26|. 
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FIG. 4. Reynolds vs Prandtl in regime 4 {Pr > 0.35, velocity fluctuations dominate but 
temperature fluctuations are negligible). The symbols is the theoretical formula. The dotted line 
is a power-law fit, with slope —0.8, mimicking the experimental fit of [26]. 
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FIG. 5. Reynolds vs Rayleigh in regime 4 {Pr > 0.35, velocity fluctuations dominate but 
temperature fluctuations are negligible). The symbols is the theoretical formula. The dotted line 
is a power-law fit, with slope 0.43, mimicking the experimental fit of [26]. 
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FIG. 6. Nusselt vs Rayleigh in regime 5 (Pr < 0.35, fluctuations dominates. The symbols 
are the power-law fits of experimental measurements by [27]. The line is the theoretical formula 
predicted by the toy model (eq. ( |20| ) with e = 0). 
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FIG. 7. Nusselt vs Rayleigh in regime 6 {Pr > 0.35, fluctuations dominates. The symbols are 
the power-law fits of the experimental measurements by pq] . The line is the theoretical formula 
predicted by the toy model (eq. ( |20| ) with e = —1/2). 
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